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An integrated computational model of a generic hypersonic vehicle has been 
developed for the purpose of determining the vehicle's performance characteristics, which 
include the Lift, Drag, Thrust, and Moment acting on the vehicle at specified altitude, flight 
condition, and vehicular configuration. The material presented herein is based on the prior 
work reported in [1], The Lift, Drag, Thrust, and Moment are developed for the body fixed 
coordinate system. These forces and moments arise from both aerodynamic and propulsive 
sources. SCRAMjet engine performance characteristics, such as fuel flow rate, can also be 
determined. 

The vehicle is assumed to be a lifting body with a single aerodynamic control surface. 
The body shape and control surface location are arbitrary and must be defined. 

The aerodynamics are calculated using either 2-Dimensional Newtonian or Modified 
Newtonian theory and approximate high-Mach-number Prandtl-Meyer expansion theory. 
Skin-friction drag has also been accounted for. The skin-friction drag coefficient is a 
function of the freestream Mach number. The data for the skin-friction drag coefficient 
values were taken from NASA Technical Memorandum 102610 [2]. 

The modeling of the vehicle's SCRAMjet engine is based on quasi 1 -Dimensional 
gasdynamics for the engine diffuser, nozzle, and the combustor with heat addition. The 
engine has three variable inputs for control. These are the engine inlet diffuser area ratio, the 
total temperature rise through the combustor due to combustion of the fuel, and the engine 
internal expansion nozzle area ratio. 

The pressure distribution over the vehicle's lower aftbody surface, which acts as an 
external nozzle, is calculated using a combination of quasi 1 -Dimensional gasdynamic theory 
and Newtonian or Modified Newtonian theory. The exhaust plume shape is determined by 
matching the pressure inside the plume, calculated from the gasdynamic equations, with the 
freestream pressure, calculated from Newtonian or Modified Newtonian theory. In this 
manner, the pressure distribution along the vehicle afterbody expansion surface is then 
determined. 

The following three sections describe the aerodynamic modeling, the engine 
modeling, and the exhaust plume analysis in more detail. Section four gives a description of 
the computer code used to perform the above calculations, and an input/output example is 
then given in section five. The computer code is available on a Macintosh floppy disk. 
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AERODYNAMIC ANALYSIS 


Compression Surface 

The non-dimensional pressure coefficient as given by Modified Newtonian theory is 
given below as [3] 


Cpcom(Y> 0 ) — 
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where 0 is the local flow deflection angle defined as 


0 = a + x 


( 2 ) 


The angles a and x are, respectively, the angle of attack and the local surface slope relative 
to the body fixed axis as shown in Fig 1. For Newtonian theory, the pressure coefficient 
value is given by Cp com (Y>M<»> 0) = 2.0 sin 2 (0). The compression surface is divided into 

sections or panels each of which has a constant value for x, say Xi. The pressure coefficient 
on the corresponding panel is then given as 

Cpi = Cp CO m(Y> M^, 0i) (3) 

where 0i = a + Xi. The force per unit width acting on the i^ 1 panel is given by 

= -Poo (1+ \ Y Moo 2 CpO dSi n, (4) 

where dsj is the area per unit width of the i l “ panel and n; is a unit vector normal to the i 1 * 1 
panel pointing outward, also shown in Fig 1. The moment per unit width acting about the 
vehicle's center of mass is given by 

Mj = x Fi (5) 

where fi is the vector from the vehicle's mass center to the centroid of the i^ 1 panel. 


Control Surface 

The force per unit width acting on the control surface is given by 


F cs = P„ 


l+iyMoo 2 C Pc 




_ /s 

^CS ^cs 


( 6 ) 


where 0 CS = a + 5, and 8 is the control surface deflection, S C s is the control surface area per 
unit width, and n cs is a unit vector normal to the control surface as shown in Fig 2. The 
moment per unit width acting about the vehicle's mass center is given by 

Me S = r cs xFcs (7) 

where r cs is a vector from the vehicle's mass center to the center of the control surface, which 
is where the aerodynamic center of the control surface is assumed to be located. 
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Expansion Surface 

The pressure coefficient as given by the approximate high Mach number Prandtl- 
Meyer expansion theory is [4] 


Cpexp(Y> CC ) — 
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where a is the angle of attack. With this, the force per unit width acting on the upper surface 
is given by 


F us = - P 


1 2 Y Cp e xp(y, M^, tx) L u n u 


(9) 


where L u is the area per unit width of the vehicle's upper surface and n u is a unit vector 
normal to the upper surface pointing outward as shown in Fig 1 . The moment per unit width 
acting about the vehicle's mass center is given by 

M us = r us x Fys (10) 


where r us is a vector from the vehicle's mass center to the center of the vehicle's upper 
surface. 


Skin-Friction Drag 

The skin-friction drag term is given as 

D sf = q 00 C do (MJS ref (11) 

where qoo is the freestream dynamic pressure. The data used for C^fM^) and S re f were 

obtained from NASA Technical Memorandum 102610 "Hypersonic Vehicle Simulation: 
Winged Cone Configuration" [2]. This drag term is in the velocity fixed coordinate system. 
In the body fixed coordinate system, the lift and drag terms are as follows 

L 0 = D s fsin(a) (12) 

D 0 = D sf cos(a) (13) 

The Total Aerodynamic Force and Moment 

The total aerodynamic force and moment per unit width is determined by adding all 

the separate contributions from the compression surface panels, the control surface, the upper 

surface, and the skin-friction drag. 

F A = lF i + F us + F cs -(D 0 i + L 0 k) (14) 

i 

M A = SM i + M us + M cs (15) 
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ENGINE ANALYSIS 


The vehicle's SCRAMjet engine consists of three main sections which are the 
diffuser, the combustor, and the internal expansion nozzle. There are four stations at which 
the engine gas flow pressure, temperature, and Mach number are calculated, as shown in Fig 
3. Station 1 is located at the engine inlet, station 2 is at the diffuser exit/combustor inlet, 
station 3 is at the combustor exit/nozzle inlet, and s.ation 4 is located at the engine exit. The 
pressure, temperature, and Mach number at each station are calculated according to the 
theory described below. The thrust produced by the engine is then calculated according to 
the following equation [1] 


Thrust = 


y[m4 2 p 4 -MliPi 


p 4 


i. 


An Ad/J 




(16) 


An Ad, 

where. Ad is the diffuser area ratio. An is the internal nozzle area ratio, and A e is the engine 
exit area per unit width. 


Engine Inlet 

At the inlet to the engine, the oncoming freestream flow is turned by the bow shock 
which, in accordance with Newtonian/Modified Newtonian theory, is parallel to the lower 
forebody surface at the inlet to the engine. The pressure, temperature, and Mach number 
behind the shock are then the inlet conditions to the vehicle's sSCRAMjet engine. The inlet 
pressure, temperature, and Mach number are given in terms of the freestream conditions and 
the engine inlet slope by the following set of equations which are consistent with 
Newtonian/Modified Newtonian theory [1] 

Pi = P„ [ 1 + y Mj Cp com (Y, M m , 0 mlet ) ] (17) 

Ti =T 00 1 + £ (Y- 1) M„ 2 sin 2 (0 inlet ) ] (18) 

M^ cos(6 in iet) 

Ml — p = ====== (19) 

Y 1 + 2 (Y-l) M ~ 2 sin 2 (0 inlet ) 

where ©inlet = a + Xiniet and Tiniet Is the inlet slope relative to the vehicle's body fixed axis as 
shown in Fig 3. 

After the flow is turned by the shock, it is assumed to run parallel with the inlet slope 
angle. This flow is then turned again, by the engine inlet, so that the flow is now parallel to 
the axis of the engine. This turning as assumed to occur without a change in the mass flow 
rate momentum; although, there is a significant turning force, which acts at the engine inlet, 
due to the change in direction of the flow velocity. This turning force per unit width can be 
calculated as [5] 
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( 20 ) 


Finiet = Y M i 2p i [(l-cose m iet)i + (sine in]et )k]-^— 

Ad An 

which can be seen from the vector diagram in Fig 4. Again, it is noted that this turning force 
is due solely to the change in the flow velocity direction as it turns into the engine. The 
moment per unit width acting about the vehicle mass center due to this turning force is given 

by _ 

Minlet = runlet X Finiet (21) 

where r m i e t is the vector from the vehicle mass center to the engine inlet. 


Engine Diffuser 

The gasdynamic equations which relate the diffuser exit conditions to the diffuser 
inlet conditions are given by [4] 
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( 22 ) 


(23) 


(24) 


The inlet conditions to the diffuser are simply the inlet conditions to the engine given by 
Eqns (17), (18), and (19). 


Engine Combustor 

The gasdynamic equations which relate the combustor exit conditions to the 
combustor inlet conditions are given by [4] 

M 3 2 l+ L (Y-l) M 3 2 M 2 2 [l+i- (Y-l) M 2 2 
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|m 2 


(25) 


(26) 


(27) 


The inlet conditions to the combustor are simply the exit conditions from the diffuser, which 
are given by Eqns (22), (23), and (24). The quantity T 0 , which denotes the total temperature 
rise across the combustor due to combustion of the fuel, is related to the mass flow rate of 
fuel by the following equation [6] 
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where m^ei is the mass flow rate of fuel to the engine, m air is the mass flow rate of air 
through the engine, r| c is the combustion efficiency, Q is the heating value of the fuel, and c p 
is the specific heat at constant pressure for the fuel/air mixture. 

In the above equation for M 3 ; i.e., Eqn (25), it is important to note that M 3 is 
essentially a function of diffuser area ratio. Ad, and total temperature rise, T 0 . For given 
values of Ad and T 0 , the exit Mach number, M 3 , must remain above M 3 = 1; i.e., the flow 
must not be choked. 


Engine Internal Expansion Nozzle 

The gasdynamic equations which relate the nozzle exit conditions to the nozzle inlet 
conditions are given by [4] 


(> 
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(29) 

(30) 

(31) 


The inlet conditions to the nozzle are simply the exit conditions from the combustor which 
are given by Eqns (25), (26), and (27). These nozzle exit conditions are also the engine exit 
conditions which are used in Eqn (16), along with the engine inlet conditions given by Eqns 
(17), (18), and (19), to determine the thrust produced by the engine. 


The Total Engine Force and Moment 

With the thrust given by Eqn (16), the thrust force per unit width is given as 

Fxh = Thrust i (32) 

The moment per unit width acting about the vehicle mass center is given by 

Mrh = (z m ] et Thrust) j (33) 

where zi n i e t is the distance from the mass center to the centerline of the engine, which is 
assumed to be parallel with the vehicle's body fixed x-axis. 

The total engine force and moment per unit width is then given by the sum of the 
separate contributions as 
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EXHAUST PLUME ANALYSIS 


w 


tj 


Plume Position and Pressure Calculation 

In Fig 5, the geometry used to perform the plume shape calculation is shown. The 
region enclosed by the vehicle's lower aftbody expansion surface and the plume surface is 
discretized into several sections. The plume shape is assumed to have a quadratic form along 
each of these sections given by 

yj = a, + b, (x - x,_i) + q (x - Xj.i) 2 (36) 

The slope of the plume shape, denoted as (3j, is then simply given by the derivative of Eqn 

(36) 


tan(pi) = bi + 2 Ci (x - x^) (37) 

The Mach number and pressure at station (i) is given in terms of the Mach number and 
pressure at station (i-1) through the gasdynamic equations for an expansion [4] 


(l + I(y-l)M 1 2 fr (l +J-(y-l)M i . 1 2 )rT 

1 2 ' = (A D )i 2 2- ' 
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(P.) inner “ (Pi-i) inner 

The area ratio for section (i) is given by 


Mu 2 


i+ky-UM^'p 
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1 +2 (Y'l)Mi 2 


(Ao)i 


_ yi(xj) - Sj 


Yi(xi-i) - Si -1 


(38) 


(39) 


(40) 


where sj is the position of the vehicle's lower aftbody surface at station (i). The pressure 
given by Eqn (39) is the pressure on the inside of the plume at station (i). The pressure on 
the outside of the plume at station (i) is given by Newtonian/Modified Newtonian theory as 


(Pi)outer — Poo 1 + Y M^ 2 Cpcomfy, M^, Gj) 


(41) 


where 0j = a + Pi. 

The plume shape at each section has three coefficients to be determined, aj, bi, and cj. 
Two of the coefficients, ai and bi, can be determined by enforcing that the plume position and 
slope at each station is continuous; i.e., 

yi(x I .i) = ai = y i .i(x i . 1 ) (42) 

<«) 

The coefficient ci is determined by imposing the condition that the inner pressure value 
determined from Eqn (39) and the outer pressure determined from Eqn (41) are equal. An 
iterative approach is used to determine the proper value of ci at each station. 
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The solution procedure begins at the engine exit, where the plume position and the 
pressure are known. The initial values for ai and bi are given by 

= 1 a (44) 


bi = tan 


"I*- 


sin 



t 7 


2 y CpcomCy, y) 


a 


(45) 


The solution at each station is then determined successively using the solution obtained at the 
previous station. 


The Total Exhaust Gas Force and Moment 

With the pressure at each station now known from the above procedure, the force per 
unit width acting on each panel of the vehicle's lower aftbody expansion surface can be 
written as 

Fi = - 2 (Pi t + ?i) ds, n, (46) 

where dsi is the area per unit width of the i l " panel and Hi is a unit vector normal to the i 1 * 1 
panel pointing outward. The moment per unit width acting about the vehicle mass center due 
to the i 1 * 1 panel is given by 

M i = f 1 xF, (47) 

where ?i is a unit vector from the mass center to the centroid of the i* panel. 

The total force and moment per unit width is then just the sum of the separate 
contributions from each of the panels 


F E x = £F, 
M E x = IH 


(48) 

(49) 



Figure 5. Plume Shape Calculation Coordinate Sysytem 






MATLAB® PROGRAM DESCRIPTION 


In this section, a program flow chart is given followed by a description of the 
subroutines called by the main program, NASPmain. The program is written in MATLAB®. 
The input file is named NASPinput and contains the information as shown in the example 
case of the next section. 

Program Flow Chart 
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Program NASPmain 

The inputs are: 

given in the input file NASPinput 
The outputs are: 

any of the following variables described below 
Variable names description: 

Cpmax - Maximum pressure coefficient value determined from Newtonian 
or Modified Newtonian Theory 
gama - Ratio of specific heats, set to 1.43 

dtr - Conversion factor for converting from degrees to radians 

rtd - Conversion factor for converting from radians to degrees 

R - Gas constant for air, set to 1686.0 [(lb ft)/(slug 0 K)] 

To - Total temperature rise across the combustor [° K] 

Ad - Engine diffuser area ratio 

Vi - Freestream velocity [ft/s] 

Machi - Freestream Mach number 

Ti - Freestream temperature [° K] 

dynpres - Dynamic pressure [lb/ft 2 ] 

Pi - Freestream pressure [lb/ft 2 ] 

arotyp - Parameter used to determine type of 'Cpmax' calculation 
arotyp = 0: Newtonian theory 
arotyp = 1: Modified Newtonian theory 
uspts - Array of two (x,z) data points defining the upper surface [ft] 

lspts - Array of (x,z) data points defining the lower forebody compression 

surface [ft] 

alpha - Angle of attack [degrees] 

Lb - Lift per unit width due to upper and lower compression 
surface[lb/ft] 

Db - Drag per unit width due to upper and lower compression 
surface [lb/ft] 

Mb - Moment per unit width about the mass center due to upper and 
lower compression surface [lb] 

delta - Control surface deflection, positive down [degrees] 

Scs - Control surface area per unit width (includes both sides) [ft] 
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xcs - x-axis location of the control surface [ft] 

zcs - z-axis location of the control surface [ft] 

Lcs - Lift per unit width due to the control surface [Ib/ft] 

Dcs - Drag per unit width due to the control surface [Ib/ft] 

Mcs -Moment per unit width about the mass center due to the control 
surface[lb] 

Laero - Total lift per unit width due to aerodynamics [lb/ft] 

Daero - Total drag per unit width due to aerodynamics [Ib/ft] 

Maero - Total moment per unit width about the mass center due to 
aerodynamics [lb/ft] 
inslp - Engine inlet angle [degrees] 
xinlt - x-axis position of the engine inlet [ft] 

zinlt - z-axis position of the engine inlet [ft] 

ctyp - parameter used in ENGINE to set the type of engine control 

ctyp = 0: calculate flow quantities based on given Ad and To 

ctyp = 1 : For the given To, Ad is calculated which will give a 
choked flow condition at the combustor exit 
ctyp = 2: For the given Ad, To is calculated which will give a 
choked flow condition at the combustor exit 
An - Engine internal nozzle area ratio 
Ae - Engine exit area per unit width [ft] 

Leng - Lift per unit width due to the engine [lb/ft] 

Teng - Thrust per unit width due to the engine [lb/ft] 

Meng - Moment per unit width about the mass center due to the engine [lb] 
Engdat - Matrix containing the engine flow quantity values 

Engdat = PI [Ib/ft 2 ] P2 [lb/ft 2 ] P3 [lb/ft 2 ] P4 [lb/ft 2 ] 

T1 [K] T2 [K] T3 [K] T4 [K] 

Ml M2 M3 M4 

exspts - Array of (x,z) data points defining the lower aftbody expansion 
surface [ft] 

Lex - Lift per unit width due to the exhaust gas [lb/ft] 

Tex - Thrust per unit width due to the exhaust gas [lb/ft] 

Mex - Moment per unit width about the mass center due to the exhaust 
gas [lb] 

plmpts - Array of (x,z) data points defining the calculated plume 
position [ft] 
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Lift 

Drag 

Thrust 

Momnt 

igraph 

scale 


- Total lift per unit width on the vehicle [lb/ft] 

- Total drag per unit width on the vehicle [lb/ft] 

- Total thrust per unit width on the vehicle [lb/ft] 

- Total moment per unit width about the mass center on the 

vehicle [lb] 

- draws the vehicle and plume shape for igraph = 1 

- sets the axis scale for the vehicle/plume diagram [ft] 


The maximum value for the compression surface pressure coefficient is first 
determined based on Newtonian or Modified Newtonian theory depending on the value of 
'arotyp'. Next, the program calls AEROBODY, which calculates the lift, drag, and moment 
values due to the vehicle's upper expansion surface, lower compression surface, and the skin- 
friction drag, and AEROCNSF, which calculates the lift, drag, and moment values due to the 
control surface, and then adds them to obtain the total aerodynamic contribution to the 
vehicle's total lift, drag, and moment. Note that the 'uspts' array contains only two data points 
which are the coordinates of the trailing edge and leading edge of the upper surface, which is 
assumed to be parallel with the body fixed x-axis. The engine contribution to the total 
vehicle lift, thrust, and moment are then calculated by calling ENGINE. Finally, PLUME is 
called which calculates the lift, thrust, and moment due to the exhaust gas flow over the 
vehicle's lower aftbody expansion surface. The plume shape is also calculated and returned 
for graphing if desired. The total vehicle lift, drag, thrust, and moment is then calculated by 
adding the separate contributions from the aerodynamics, engine, and exhaust gas flow. 

Subroutine MODCP 

The inputs are: 

M - Freestream Mach number 


I 


The outputs are: 

cpvalue - The maximum pressure coefficient value calculated from 
Modified Newtonian theory 

This subroutine calculates the maximum pressure coefficient value based on Modified 
Newtonian theory which is the coefficient on the sin 2 (0) term in Eqn (1). 
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Subroutine AEROBODY 

The inputs are: 

uspts - Array of (x,z) data points defining the upper surface [ft] 

Ispts - Array of (x,z) data points defining the lower forebody compression 

surface [ft] 

Mi - Freestream Mach number 

Pi - Freestream pressure [lb/ft 2 ] 

qi - Freestream dynamic pressure [lb/ft 2 ] 

alpha - Angle of attack [degrees] 

The outputs are: 

L - Lift per unit width due to upper and lower compression 

surface[lb/ft] 

D - Drag per unit width due to upper and lower compression 

surface [Ib/ft] 

M - Moment per unit width about the mass center due to upper and 
lower compression surface [lb] 

The subroutine first determines if the upper surface is an expansion surface (a > 0) or 
a compression surface (a < 0) and determines the upper surface pressure accordingly. The 
two points defining the upper surface are then reordered, through a call to REORDR, so that 
the GTNORML subroutine returns the correct direction for the upper surface unit normal 
vector. The upper surface force and moment is then calculated according to Eqns (9) and 
(10). The same procedure is followed for the lower compression surface panels. For each 
panel, the unit normal vector is determined from which the panel force and moment are 
calculated according to Eqns (4) and (5). The separate panel forces and moments are then 
summed to obtain the total force and moment due the compression surface. Next, the skin- 
friction drag coefficient is obtained for the particular freestream Mach number input. The 
skin-friction drag force is then calculated and transformed from the velocity fixed coordinate 
system to the body fixed coordinate system using Eqns (1 1), (12), and (13). The total lift and 
drag is then determined from the components of the total force vector. 

Subroutine AEROCNSF 

The inputs are: 

delta - Control surface deflection, positive down [degrees] 

Scs - Control surface area per unit width (includes both sides) [ft] 


m 
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- x-axis location of the control surface [ft] 

- z-axis location of the control surface [ft] 

- Freestream Mach number 

- Freestream pressure [lb/ft 2 ] 

- Dynamic pressure [lb/ft^] 

- Angle of attack [degrees] 

- Lift per unit width due to the control surface [Ib/ft] 

- Drag per unit width due to the control surface [lb/ft] 

-Moment per unit width about the mass center due to the control 

surface [lb] 

This subroutine calculates the lift, drag, and moment about the vehicle mass center 
due to the control surface. The force acting on the control surface is first calculated in the 
velocity fixed coordinate system. This force is then transformed to the body fixed coordinate 
system depending on which side (upper or lower) of the control surface is exposed to the 
oncoming freestream flow. The lift and drag are then determined from the components of the 
force vector. The moment about the vehicle mass center is calculated according to Eqn (7). 

Subroutine ENGINE 

The inputs are: 

alpha - Angle of attack [degrees] 
inslp - Engine inlet angle [degrees] 
xinlt - x-axis position of the engine inlet [ft] 

zinlt - z-axis position of the engine inlet [ft] 

Pi - Freestream pressure [lb/ft 2 ] 

Ti - Freestream temperature [* K] 

Mi - Freestream Mach number 

ctyp - Parameter used in ENGINE to set the type of engine control 

ctyp = 0: calculate flow quantities based on given Ad and To 
ctyp = 1 : For the given To, Ad is calculated which will give a 
choked flow condition at the combustor exit 
ctyp = 2: For the given Ad, To is calculated which will give a 
choked flow condition at the combustor exit 
An - Engine nozzle area ratio 


xcs 

zcs 

Mi 

Pi 

qi 

alpha 


The outputs are: 
L 

D 

M 


18 



Ae - Engine exit area per unit width [ft] 

The outputs are: 

L - Lift per unit width due to the engine [Ib/ft] 

T - Thrust per unit width due to the engine [Ib/ft] 


My 

- Moment per unit width about the mass center due to the 

engine [lb] 

data 

- Matrix containing the engine flow quantity values 



data = PI [lb/ft 2 ] 

P2 [lb/ft 2 ] 

P3 [Ib/ft 2 ] 

P4 [lb/ft 2 ] 


T1 [K] 

T2[K] 

T3 [K] 

T4 [K] 


Ml 

M2 

M3 

M4 


PI - Pressure at engine/diffuser inlet 
T1 - Temperature at engine/diffuser inlet 
Ml - Mach no. at engine/diffuser inlet 
P2 - Pressure at diffuser exit/combustor inlet 
T2 - Temperature at diffuser exit/combustor inlet 
M2 - Mach no. at diffuser exit/combustor inlet 
P3 - Pressure at combustor exit/nozzle inlet 
T3 - Temperature at combustor exit/nozzle inlet 
M3 - Mach no. at combustor exit/nozzle inlet 
P4 - Pressure at nozzle/engine exit 
T4 - Temperature at nozzle/engine exit 
M4 - Mach no. at nozzle/engine exit 

The inlet conditions to the engine are first calculated using Eqns (17), (18), and (19). 
Next, depending on the type of engine control; i.e. 'ctyp' value, the flow conditions across the 
diffuser, combustor, and internal nozzle are calculated through calls to the subroutines 
DIFFUSER, COMBUSTOR, NOZZLE, GTAD, and GTTO. The subroutines GTAD and 
GTTO are called to determine the diffuser area ratio. Ad, or total temperature rise required, 
T 0 , to give a choked flow condition; i.e., M3 = 1.1, at the combustor exit. The values are 
then stored in the matrix 'data'. The engine thrust and turning force are then determined using 
Eqns (16) and (20) respectivley. The total moment about the vehicle's mass center is then 
calculated from Eqn (35). 
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Subroutine PLUME 

The inputs are: 
exspts 


Pi 

Mi 

Engdat 


Ae 

alpha 


- Array of (x,z) data points defining the lower aftbody expansion 
surface [ft] 


- Freestream pressure [lb/ft 2 ] 

- Freestream Mach number 

- Matrix containing the engine flow quantity values 

Engdat = P 1 [lb/ft 2 ] P2 [lb/ft 2 ] P3 [lb/ft 2 ] 

T1 [K] T2 [K] T3 [K] 

Ml M2 M3 

- Engine exit area per unit width [ft] 

- Angle of attack [degrees] 


P4 [lb/ft 2 ] 
T4 [K] 
M4 


The outputs are: 

L - Lift per unit width due to the exhaust gas [lb/ft] 

T - Thrust per unit width due to the exhaust gas [lb/ft] 

My - Moment per unit width about the mass center due to the exhaust 
gas [lb] 

plmpts - Array of (x,z) data points defining the calculated plume 
position [ft] 



The lower aftbody surface points are first transformed to the local plume calculation 
coordinate system as shown in Fig 5 and the initial values for the coefficients 'a' and 'b' are 
set according to Eqns (44) and (45). 

The iterative process for determing the plume position begins by initialy setting the 
coefficient 'c' to c = 0. The section area ratio is then determined form Eqn (40) and an inner 
iteration method is used to solve Eqn (38) for the current station Mach number. The inner 
pressure is then determined from Eqn (39). The slope of the plume is calculated from Eqn 
(37) and used in Eqn (41) to determine the outer pressure due to aerodynamics. A Newton 
Method is then used to obtain an updated value for 'c' based on the difference in inner and 
outer pressure value. The above procedure is repeated until the difference in inner and outer 
pressure values is within a specified tolerance. The (x,z) data points defining the position of 
the plume are then transformed back to the body fixed coordinate system and stored in the 
array 'plmpts'. 
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The force on each lower surface panel is determined according to Eqn (46) and the 
total lift and thrust is determined from the total force vector components. The total moment 
about the vehicle mass center is determined according to Eqns (47) and (49). 

Subroutine DIFFUSER 

The inputs are: 

PI - Diffuser inlet pressure [lb/ft 2 ] 

T1 - Diffuser inlet temperature [° K] 

Ml - Diffuser inlet Mach number 

The outputs are: 

P2 - Diffuser exit pressure [lb/ft 2 ] 

T2 - Diffuser exit temperature [° K] 

M2 - Diffuser exit Mach number 

The subroutine simply solves Eqns (22), (23), and (24) for the diffuser exit conditions 
for a specified diffuser area ratio. An iterative method is used to solve Eqn (22) for 'M2'. 

Subroutine COMBUSTOR 

The inputs are: 

PI - Combustor inlet pressure [lb/ft 2 ] 

T1 - Combustor inlet temperature [° K] 

Ml - Combustor inlet Mach number 

The outputs are: 

P2 - Combustor exit pressure [lb/ft 2 ] 

T2 - Combustor exit temperature [° K] 

M2 - Combustor exit Mach number 

The subroutine simply solves Eqns (25), (26), and (27) for the combustor exit 
conditions for a specified total temperature rise across the combustor. The subroutine checks 
for choked flow conditions. If the inlet conditions and total temperature rise value 
combination yields a subsonic exit Mach number, then an error message is displayed. An 
iterative method is used to solve Eqn (25) for 'M2'. 
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Subroutine NOZZLE 


The inputs are: 


PI 

- Nozzle inlet pressure [lb/ft 2 ] 

T1 

- Nozzle inlet temperature [° K] 

Ml 

- Nozzle inlet Mach number 

An 

- Engine nozzle area ratio 


The outputs are: 

P2 - Nozzle exit pressure [lb/ft 2 ] 

T2 - Nozzle exit temperature [* K] 

M2 - Nozzle exit Mach number 

The subroutine simply solves Eqns (29), (30), and (31) for the internal nozzle exit 
conditions. An iterative method is used to solve Eqn (29) for 'M2'. 

Subroutine GT AD 

The inputs are: 

PI - Diffuser inlet pressure [lb/ft 2 ] 

T1 - Diffuser inlet temperature [° K] 

Ml - Diffuser inlet Mach number 

The outputs are: 

P2 - Diffuser exit/Combustor inlet pressure [lb/ft 2 ] 

T2 - Diffuser exit/Combustor inlet temperature [° K] 

M2 - Diffuser exit/Combustor inlet Mach number 

P3 - Combustor exit pressure [lb/ft 2 ] 

T3 - Combustor exit temperature [° K] 

M3 - Combustor exit Mach number 


This subroutine determines the diffuser area ratio required for combustor choked flow 
exit conditions given engine/diffuser inlet conditions and the total temperature rise through 
the combustor. Eqn (24) is first substituted into Eqn (25) for T2‘, then the resulting equation 
is solved for 'M2' which gives a choked flow combustor exit Mach number M3 =1.1. With 
'M2' known, Eqn (22) is used to determine 'Ad', the diffuser area ratio. The diffuser exit 
pressure and temperature are determined from Eqns (23) and (24) and the combustor exit 
pressure and temperature are determined from Eqns (26) and (27). 
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Subroutine GTTO 







The inputs are: 

PI - Combustor inlet pressure [lb/ft 2 ] 

T1 - Combustor inlet temperature [° K] 

Ml - Combustor inlet Mach number 

The outputs are: 

P2 - Combustor exit pressure [lb/ft 2 ] 

T2 - Combustor exit temperature [° K] 

M2 - Combustor exit Mach number 

This subroutine determines the total temperature rise through the combustor required 
for combustor choked flow exit conditions given the combustor inlet conditions. Eqn (25) is 
solved for 'To', the total temperature rise through the combustor. The choked flow exit Mach 
number is set at M2 = 1.1. Eqns (26) and (27) are used to determine the exit pressure and 
temperature. 

Subroutine REORDR 

The inputs are: 

xypts -An array containing (x,z) data points to be reordered 
sgn - A parameter specifying the type of reordering to be performed 
sgn > 0 points ordered such that the x coordinates are in 
ascending order 

sgn < 0 points ordered such that the x coordinates are in 
descending order 

The outputs are: 

xyptso -An array containing the reordered (x,z) data points 

This subroutine simple orders the (x,z) data points such that the x coordinate is in 
ascending or descending order, depending on the sign of 'sgn'. 

Subroutine GTNORML 

The inputs are: 

x2 - x coordinate of data point (x,z)2 

z2 - z coordinate of data point (x,z)2 


■ 
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xl 

zl 


- x coordinate of data point (x,z) i 

- z coordinate of data point (x,z)i 


The outputs are: 

r - vector from the coordinate system origin to the center of a line 

connecting points (x,z)i and (x,z )2 

n - unit vector normal to the line pointing away from the origin, 

ds - length of the line 

This subroutine first determines the components of the vector 'r* and the length of the 
line 'ds', from the input data points. The unit normal is then calculated according to the 
following set of equations 

n x t = - 1 1 1 
~n • t = 0 

* 

where t is the vector from (x,z)i to (x,z) 2 . 



INPUT/OUTPUT EXAMPLE 


»_ j 

m 

m 


m 

H 

y 


t3 

m 


In this section an example case is given. The input file, NASPinput, is shown below. 
The output forces and moments are given on pages 24 and 25. 

Sample NASPinput File 

% Flight Conditions at 100,000 ft altitude 
Machi=8; 

Pi=44.95; 

Ti=396.36; 

% 

% Vehicle Configuration 

alpha=2; 

delta=5; 

Scs=22.5; 

xcs=-35; 

zcs=-10; 

% 

% Engine Configuration 
inslp= 18.4349; 

Ae=5; 

xinlt=5; 

zinlt=10; 

Ad=0.1482; 

To=2000.0; 

An=6.35; 

% 

% Newtonian (arotyp=0) or Modified Newtonian (arotyp=l) 
arotyp=0; 

% 

% Type of engine control 
% ctyp=0 - straight calculation 

% ctyp=l - determine Ad to choke the flow for given To 
% ctyp=2 - determine To to choke the flow for given Ad 
ctyp=0; 
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% 


% Body Geometry 
uspts=[-35, -10; 
65, -10]; 


lspts=[65.0000, -10.0000; 


59.0000, 

-8.0000; 

53.0000, 

-6.0000; 

47.0000, 

-4.0000; 

41.0000, 

-2.0000; 

35.0000, 

-0.0000; 

29.0000, 

2.0000; 

23.0000, 

4.0000; 

17.0000, 

6.0000; 

11.0000, 

8.0000; 

5.0000, 

10.0000; 

-5.0000, 

10.0000]; 


% 

exspts=[-5.0000, 10.0000; 


-8.0000, 

6.6761; 

- 11 . 0000 , 

3.6464; 

-14.0000, 

0.9109; 

-17.0000, 

-1.5304; 

-20.0000, 

-3.6775; 

-23.0000, 

-5.5304; 

-26.0000, 

-7.0891; 

-29.0000, 

-8.3536; 

-32.0000, 

-9.3239; 

-35.0000, 

-10.0000]; 


igraph=0; 

scale=65; 
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Sample NASP output 


Lb = 3.2443e+04 
Db = 7.4344e+04 
Mb = 1.1016e+06 

Lcs = 2.3770e+03 
Dcs = 207.9627 
Mcs = -8.1 1 16e+04 

Laero = 3.4820e+04 
Daero = 7.4552e+04 
Maero = 1.0205e+06 

Leng = -2.7557e+04 
Teng = 1.4958e+04 
Meng = 1.1791e+04 

Lex = 5.9956e+03 
Tex = 5.2346e+03 
Mex = -6.0108e+04 

Lift = 1.3258e+04 
Drag = 7.4552e+04 
Thrust = 2.0192e+04 
Momnt = 9.7221e+05 
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